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The gauge group being centreless, G2 gauge theory is a good laboratory for studying the role of the 
centre of the group for colour confinement in Yang-Mills gauge theories. In this paper, we inves- 
tigate G2 pure gauge theory at finite temperature on the lattice. By studying the finite size scaling 
of the plaquette, the Polyakov loop and their susceptibilities, we show that a deconfinement phase 
transition takes place. The analysis of the pseudocritical exponents give strong evidence of the 
deconfinement transition being first order. Implications of our findings for scenarios of colour 
confinement are discussed. 
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1. Introduction and Motivation 

Quark confinement is one of the oldest issues in non perturbative QCD. The fundamental prob- 
lem is the characterization of the degrees of freedom responsible of confinement. Many choices 
are investigated in literature, among them center vortices and abelian monopoles (both proposed by 
't Hooft) are the most popular. The connection between centre symmetry breaking and deconfine- 
ment transition in SU (N) theories, and the Svetitsky-Yaffe conjecture [p relating the universality 
class of the deconfinement transition of a gauge theory to that of one with its centre as a spin sistem 
in one lower dimension, seem to link the presence of a non trivial centre to confinement. It's there- 
fore an interesting question whether centerless groups present or not a confmement-deconfinement 
phase transition. One example in this direction is the 50(3) = SU (2)/Z2 that has been extensively 
studied on the lattice (see for example [Q, |3|, Another choice is the exceptional group G2 
which moreover is without 't Hooft center vortices and is the simplest one with this property, so we 
focused on it for our investigations. 

Another point of view concerning the confinement mechanism is the so called dual super- 
conductor picture, the subject of the project. The spontaneous breaking of the electromagnetic 
symmetry to a discrete group in ordinary superconductors, due to condensation of Cooper pairs, 
leads to the formation of the so called Abrikosov flux tubes connecting magnetic charges. The dual 
mechanism should connect quark pairs by strings of (chromo)electric flux tubes due to condensa- 
tion of magnetic objects. According to this Dual Superconductor Picture (DSP), the QCD vacuum, 
is a condensate of magnetically charged fields (monopoles) confining (chromo)electrically charged 
particles. The G2 group admits stable monopole solutions, so in principle there could be a phase 
where these condense giving rise to confinement even without the presence of a centre. 

We started our investigation by analyzing in detail the fhermodynamical properties of the ex- 
ceptional group G2 and we report here our findings. The first and essential issue is to check whether 
the transition found in [||] is a real one or a crossover, for example. Were this the case, the physics 
of deconfinement in G2 would be noticeably different from that of SU (N) gauge theories, and this 
would cast serious doubts about what can we learn from G2 for confinement in more physical gauge 
theories. 

We perform a Finite Size Scaling (FSS) analysis of the plaquette and of the Polyakov loop 
susceptibilities. We test the educated guess of a first order phase transition and we can show that a 
real transition takes place so the confined and deconfined phase are really different. In section [O] 
we discuss some basic properties of the G2 group. Then, in section || we state the results of our 
lattice simulations. We draw the conclusions in the last section. 

1.1 G2 group 

We are now going to state some basic facts about the Lie Group G2. It can be naturally 
constructed as a subgroup of the real group S0{1) which has 21 generators and rank 3. To the 
usual properties of S0{1) matrices 



detQ = 1 



or 1 =o. T 



(l.i) 



we have in addition another constraint 



Tabc — Tdef^da^eb^fc 



(1.2) 
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where T a j, c is a totally antisymmetric tensor whose nonzero elements are (using the octonion basis 
given by [@]) 

^123 = ^176 = ^145 = T 257 = T 24 (, = T 3 A1 = T 3(t5 = 1. (1.3) 



Equations 1.2 are 7 relations reducing the numbers of generators to 14. The fundamental represen- 
tation of G 2 is 7 dimensional and by using the algebra representation of [||] we can clearly identify 
an SU (3) subgroup and several SU (2) subgroups, with 6 of them we can cover the whole group. 
The first three SU (2) subgroups are the 4x4 real representations of the group while the remaining 
three are extremely difficult to simulate with standard techniques. See next section for details on 
simulations. 

The Lie group G 2 has rank 2 as SU (3), this implies that its Cartan subgroup, that is the maximal 
residual abelian subgroup after an abelian gauge fixing, is U(l) X C/(l). Stable monopole solutions 
are classified according to the homotopy group 1 : 

K 2 (G 2 /U(l) 2 ) = 7Ti(C/(l) x 1/(1)) = ZxZ (1.4) 

i.e. we have two distinct species of monopoles like in SU (3). 

Another interesting homotopy group shows that center vortices are absent in the theory: 

k 1 (G 2 /^(G 2 )) = k 1 {G 2 )=0 (1.5) 

while for SU (3) for example 

K l (SU(3)/Z 3 )=Z 3 . (1.6) 

So G 2 is a good playground to study the dual superconductor picture in a theory without center 
vortices, thus isolating monopole contribution in confinement. 



2. Simulation and Results 

Much work has been devoted to simulations of G 2 Yang-Mills theory by the Bern group || HJ. 
Here we are going to investigate in some detail the thermodynamical properties of the system to 
demonstrate that a deconfmement phase transition is indeed taking place. 

To simulate the gauge theory 

JSf=r^TrF MV F MV (2.1) 

discretized with the Wilson action, we used a simple Cabibbo-Marinari update (heat-bath + overre- 
laxation in a tunable ratio) for the first three SU (2) subgroups (4x4 representations) spanning the 
SU (3) C G 2 . To guarantee the covering of the whole gauge group we make a completely random 
gauge transformation every n updates (tipically 2). The observables measured are the standard 
plaquette and the Polyakov Loop. To study the thermodynamical properties we simulated several 
lattices of spatial dimension N s = 12, 16, 18, 20 and TV, = 6 (N t = 4 only for the smallest lattice). At 
the transition we needed histories of the order of 10 5 updates. The code uses only real algebra and 
run on an Opteron farm here in the computer facilities of the Physics Department of the University 
of Pisa. 

'The first equality follows from %\ (G2) = 0. See for example [M]. 
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Figure 1: Left: plaquette susceptibility for various spatial volumes and N t . The huge peak at j3 = 1 .35 is the 
bulk transition described in the text. In the inset is shown the real phase transition point. Right: comparison 
of plaquette for non zero and zero temperature simulations. The integral of the difference of the curves gives 
the free energy density which is different from zero only at the physical transition. 



2.1 Thermodynamics 

The measured observables are the plaquette and the Polyakov loop, and their susceptibilities. 
The lattice specific heat has a strong peak at j8 ~ 1.35. This does not signal a real phase transition. 
It is present for every volume and different ./V, (a physical transition should move to the weak 
coupling region with N t ) and also at zero temperature simulations and can be ascribed to lattice 
artifacts (see Fig. [[]). This big peak completely overshadows the real physical transition that can 
be seen as a little peak in the weak coupling region at j8 ~ 1.395 for N t = 6 and several spatial 
volumes. We need a subtraction of this background to study the finite size scaling of the plaquette 
susceptibility. We estimate the background by mean of zero temperature simulations at different 
spatial volumes (to keep under control possible systematic errors). The scaling of the susceptibility 
assuming a first order phase transition is shown in Fig. ||. The specific heat scales nicely as it is 
confirmed by the linear fit of the behaviour of the plaquette susceptibility peak with the volume 




Figure 2: Left: scaling of the plaquette susceptibility with a first order ttansition hypotesis after the subtrac- 
tion of the unphysical background. Right: Linear fit of the peak heights, y = a-x + b, a = 0.00079(14) ■ 10~ 3 , 
/> = 0.98(62) -10- 3 ,;^ = 1 .35. 
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0,2 

Pol. Loop 



Figure 3: Normalized densities of the Polyakov Loop in a semilog plot for j3 varying in the range from 1 .35, 
the critical coupling of the bulk transition "Pbulk \ t0 1-401, in the deconfined phase (data from the 6 x 16 3 
lattice for the upper graph and from 6 x 20 3 for the other - same scales and limits for both axes are used 
for better comparison). As an aside we notice that far in the confined phase, j3 c < 1.395, the Polyakov loop 
is zero within errors and this feature can not be explained on the ground of any manifest symmetry of the 
system. 



The Polyakov loop is insensitive to the bulk transition. We studied its FSS even if it is not a 
proper order parameter. It has anyway a divergence at the transition, signal of a non trivial overlap 
with the real (unknown) order parameter, so that we can safely extract the critical exponents. It 
develops an evident double peak structure (see Fig. ||). From the behaviour of this observable one 
can guess a first order transition. A scaling analysis shows clearly the expected behaviour for large 
volumes, confirming the results of the specific heat. 

2.2 Conclusions 

We simulated the G2 Yang-Mills theory at finite temperature. In order to investigate the prop- 
erties of the vacuum by mean of a magnetic charged operator we need to confirm that a real 
confinement-deconfinement transition is taking place. We showed that the specific heat and the 




Figure 4: Left: Scaling of the Polyakov loop assuming first order. For the smallest lattices corrections to the 
scaling are evident; f3 c = 1.395. Right: Scaling of the peak of the Polyakov Loop susceptibility. The solid 
line is a linear fit to the data. 
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Polyakov Loop observables agree with the hypotesis of a first order transition after subtraction of 
a background due to lattice artifacts. The centre of the group seems an unessential property for 
the presence of a deconfinement transition, at least in this case. Future developements concern 
the measure of the magnetic order parameter proposed by the Pisa group to study the topological 
structure of the vacuum of this centreless theory. 

The work of CP. has been supported in part by contract DE-AC02-98CH1-886 with the U.S. 
Department of Energy. 
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